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Abstract 

The notion of microscopic state of the system at a given moment of time as a 
point in the phase space as well as a notion of trajectory is widely used in classical 
mechanics. However, it does not have an immediate physical meaning, since arbitrary 
real numbers are unobservable. This notion leads to the known paradoxes, such as the 
irreversibility problem. A "functional" formulation of classical mechanics is suggested. 
The physical meaning is attached in this formulation not to an individual trajectory 
but only to a "beam" of trajectories, or the distribution function on phase space. The 
fundamental equation of the microscopic dynamics in the functional approach is not 
the Newton equation but the Liouville equation for the distribution function of the 
single particle. The Newton equation in this approach appears as an approximate 
equation describing the dynamics of the average values and there are corrections to the 
Newton trajectories. We give a construction of probability density function starting 
from the directly observable quantities, i.e., the results of measurements, which are 
rational numbers. 

1 Introduction 



The conventional widely used concept of the microscopic state of the system in classical 
Newtonian mechanics [1] at some moment of time as the point in phase space, as well as 
the notion of trajectory and the microscopic equations of motion have no direct physical 
meaning, since arbitrary real numbers are unobservable. Observable physical quantities are 
only presented by rational numbers [2-4], see also the discussion of concepts of space and 
time in [2-8]. 

In [9] it was suggested a "functional" formulation of classical mechanics. The fundamental 
equation of the microscopic dynamics in the functional approach is not the Newton equation, 
but the Liouville equation for the distribution function of a single particle. The Newton 
equation in this approach appears as an approximate equation describing the dynamics of 
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the average values, and there are corrections to the Newton trajectories. The functional 
formulation of classical mechanics gives also an approach to the solution of the irreversibility 
problem. 

In this note we give a construction of the probability density function starting from the 
directly observable quantities, i.e., the results of measurements, which are rational numbers. 

2 States and observables in 

functional classical mechanics 

Usually in classical mechanics the motion of a point body is described by a trajectory in the 
phase space, i.e., the values of the position and momentum as functions of time, which are 
solutions of the equations of Newton or Hamilton. 

However, this mathematical model is an idealization of the physical process, rather far 
separated from reality. Every physical body has the spatial dimensions, such a mathematical 
point gives only an approximate description of the physical body. The mathematical notion 
of a trajectory does not have a direct physical meaning, since it uses arbitrary real numbers, 
i.e., infinite decimal expansions, while the observation is only possible, in the best case, in 
rational numbers, and even among them only with some error. Therefore, in the "functional" 
approach to classical mechanics, we are not starting from the Newton equation, but with 
the Liouville equation. 

Consider the motion of a classical particle along a straight line in the potential field. The 
general case of many particles in the 3-dimensional space is discussed below. Let (q,p) be 
coordinates on the plane R 2 (phase space), t G R is time. The state of a classical particle at 
time t will be described by a function p = p(q,p,t), it is the density of the probability that 
the particle at time t has the position q and momentum p. 

The description of a mechanical system with the help of probability distribution function 
P = p(Qi V) t) does not necessarily mean that we are dealing with a set of identically prepared 
ensemble of particles. Usually in probability theory one considers an ensemble of events and a 
sample space. But we can use the description with the function p = p(q,p, i) also for individ- 
ual bodies, such as planets in astronomy (the phase space in this case the 6-dimensional). In 
this case one can think on the "ensemble" of different astronomers which observe the planet, 
or on the "ensemble" of different models of behaviour of a given object for one "intelligent" 
observer. Actually, it is implicitly always dealt with the function p = p(q,p,t) which takes 
into account the inherent uncertainty in the position and momentum of the body. 

The specific type of function p depends on the method of the preparation of the state of 
the classical particle at the initial time and the type of potential field. When p = p(q,p,t) 
has sharp peaks at g = go and p = po, we say that the particle has the approximate values 
of the position and momentum q and po- 

Emphasize that the exact determination of the position and momentum can not be done 
not only in quantum mechanics, where there is the Heisenberg uncertainty relation, but also 
in classical mechanics. Always there are some errors in setting the position and momentum. 
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There are classical uncertainty relations: 



AqAp > 0, 

i.e., the uncertainty (errors of observation) in the determination of the position and mo- 
mentum is always positive (nonzero). The concept of arbitrary real numbers, given by the 
infinite decimal series, is a mathematical idealization, such numbers cannot be measured in 
the experiment. 

Therefore, in the functional approach to classical mechanics the concept of precise trajec- 
tory of a particle is absent, the fundamental concept is a distribution function p = p(q,p,t), 
and 5-function as a distribution function is not allowed. 

We assume that the continuously differentiable and integrable function p = p(q, p, t) 
satisfies the conditions: 

p>0, / p(q,p,t)dqdp=l, t ER. (1) 

If / = f(q,p) is a function on the phase space, the average value of / at time t is given 
by the integral 

7(f) = / f(q,p)p(q,P,t)dqdp. (2) 



In a sense, we are dealing with a random process £ (t) with values in the phase space. Motion 
of a point body along a straight line in the potential field will be described by the equation 

dp _ P dp | dV(g) dp 
dt m dq dq dp 

Here V(q) is a potential field and m > is the mass of the body. 

Equation (j3J) looks like the Liouville equation, which is used in statistical physics to 
describe a gas of particles, but here we use it to describe a single particle. 

The characteristics equations for fl3]) are Hamilton's equations 

. _ d_H_ . _dH_ 

^ dp ' ^ dq , 

where the Hamiltonian is 

H=^ + V(q). (5) 

Emphasize again that the Hamilton equations (jlj) in the current functional approach to me- 
chanics do not describe directly the motion of particles, but they are only the characteristics 
equations for the Liouville equation Q. 

If the distribution po{q,p) for t = is known, we can consider the Cauchy problem for 
the equation (JSJ): 

p\t=o = Po(q,p) ■ (6) 
Consider the case when the initial distribution has the Gaussian form: 

1 _(q-q ) 2 (p-pq) 2 

Po{q,P) = ^ e ~^^e . (7) 
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At sufficiently small values of the parameters a > and b > the particle has the position 
and momentum close to the go and po- For this distribution the average values of the position 
and momentum are 



Q = J qpo(q,p)dqdp = q , P= / ppo(q,p)dqdp = p 

and the dispersion 

lo . n ^ 1 



Ag 2 = (g - g) 2 = -a 2 , Ap 2 = (p - pf = -b 2 . (9) 



For the free motion (V = 0) we get 



q(t) = q + -t, p(t)=p , (10) 
m 



and the dispersion increases with time: 

Ag 2 (t) = i(a 2 + ^)- (11) 

Even if the particle was arbitrarily well localized (a 2 is arbitrarily small) at t = 0, then 
at sufficiently large times t the localization of the particle becomes meaningless, there is a 
delocalization of the particle which accounts for irreversibility 

Corrections to the Newton's trajectories for the nonlinear coupling are computed in [9]. 



3 Probability density function and rational numbers 

The probability density function is real-valued. However, our initial point was that the 
real numbers are unobservable. Does the use of real-valued probability density function 
as a fundamental notion of mechanics contradict to our initial thesis? In this section we 
construct the probability density function (formula (j!4p ) starting from the directly observable 
quantities using the methods of mathematical statistics. An important point is that the 
probability might be a real number and this is admissible, since the probabilities are not 
directly observable. 

The directly observable quantities are the results of measurements and they are rational 
numbers. Consider a measurement of an observable X. For simplicity we consider a one- 
dimensional observable, the generalization is simple. Every measurement device has an error, 
which must be taken into account. Roughly speaking, the measurement errors can be divided 
into two types: systematic and random errors [10]. If we perform repeated measurements, 
the systematic error does not change and random error changes randomly (we do not consider 
the part of the error that changes regularly, because it can be excluded by the statistical 
methods). It is natural to model the random error by a random variable. There is no 
theory of systematic error ("In fact, the only theory of systematic errors is that they must 
be identified and reduced until they are much less than the required precision... However, 
this goal is often not attainable" [10]). Therefore, although the systematic error is constant, 
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our ignorance of the systematic error is also modeled by a random variable, because of the 
absence of another theory. 

Thus, the result of a measurement is a random variable X, it is rational-valued. Moreover, 
since the precision (sensitivity) of every instrument is finite, X takes values not on the whole 
field of rational numbers Q, but rather on the lattice, X e -Z, where the rational fraction 
| is the measuring sensitivity of the instrument and Z are integers. So, the probabilities 
p m = Pi[X = |m], m G Z, are defined. The probabilities p m are real and this is admissible, 
since the probabilities are not directly observable. They can be considered as a limit of 
relative frequencies: 



p m = lim — ^ in probability, i.e., lim Pr 



Pm 

n 



0, 



where n is the number of experiments and n rn is the number of experiments where the 
realization of X is equal to |m. This is the law of large numbers [11]. The fact that 
probabilities are real numbers, actually, is not surprising, since the limit of rational sequences 
is not necessarily a rational number. 

Now consider the dynamics. Let us measure the observable X once again at some moment 
of time t > 0. We want to predict the probabilities of the results of this measurement on the 
condition that we know the result of the measurement at time t — 0. If we describe a state 
as a sum of delta functions and solve the Liouville equation with such initial conditions (this 
is equivalent to Newton's equation), we will get incorrect predictions. For example, consider 
the free motion on the real line. Assume that at time t = we obtained that momentum is 
equal to zero with the precision allowed by our instrument. Then we can conclude that at any 
time t > the particle still will be in its initial position. But in general this not true, since 
the momentum can be very small (smaller than our measuring sensitivity), but not zero. In 
this case, if t is large enough, the position of the particle can be changed considerably. 

Thus, in order to take the growth of the error with the time into account, we must 
consider the states as continuous distributions. Let us assign^ some continuous real- valued 
random variable X to our discrete random variable X. Let X be distributed according to 
some probability density function p*(x) which satisfies the condition 

Pm= p*(x)dx. (12) 

q v 2 7 

We assume that X is normally distributed: 

p ' (x)= vkt e ~^' (13) 

where x* is a mean value (which can, but not necessarily, be referred as a "true" value of 
the observable X) and a 2 is a dispersion. 

2 2 2 

° = °syst + °rand> 

where of yst and of and are the summands that correspond to the systematic and random 
error accordingly. In fact, the further discussion does not depend critically on the form of 
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distribution. We made an assumption about normal distribution for simplicity, but there 
are also some physical and mathematical reasons to choose this distributions among others. 

Again, the notion of real-valued probability density function p* does not contradict to 
the thesis that real values are unobservable, because the probability density function is 
not an observable. This is an abstract, theoretical object, which is useful, because we 
can approximate the relative frequencies using the notion of real-valued probability density 
function: 

n Ja 

Here k{X 6 [a, b]} is the number of experiments where the realization of X belongs to [a, b] 
(for example, a, b e Q) and n is the general number of experiments (it is assumed that n is 
large). 

Usually we do not now the expectation value x* and the dispersion of random error of and 
(and hence, we do not know the probability density function p*, we only assume that it has 
the form (|T3|) with unknown parameters), but rather we have to estimate them using the 
methods of mathematical statistics. The dispersion of systematic error <jg yst is assumed to 

be known from the measuring instrument certificate. Let j' 1 ', . . . ,X^ n ' be n copies of X, 
i.e., independent and identically distributed (according to the probability density function 
p*) random variables (the results of n measurements). Then the following formulas are used 
to estimate the expectation and dispersion of the random error: 

n ^— ' n — \ 

i=l i=l 

It is well-known that random variable X is normally distributed with the same expecta- 
tion x* as every of X^ 1 ', . . . ,X^ n \ The dispersion of the random error in the estimation of 
X is reduced by n times and is equal to of and /n. Hence, its estimation is S^ and /n [11]. The 
dispersion of the systematic error does not depend on the number of measurements and still 
equals to a^ yst (this is not a rigorous conclusion, since there is no theory of systematic error 
and the use of the formalism of random variables is not very correct for this). Therefore, 
the dispersion of the general error of the estimation of x* is 

C2 

q2 _ °rand , 2 
° — r O t- 

n J 

Now we construct the probability density function: 

p - (x) = vhf' e ^ ■ (14) 

If n is large (in practice, n > 30 is enough), then p n (x)Ax has the meaning of the probability 
for the mean value x* to belong to the interval (x — Ax, x + Ax). 

Remark. The last assertion can be understood by physicists, but is not completely correct 
from the mathematical point of view. Since x* is not a random variable, the probability for 
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it to belong to the definite interval is whether zero or one. More rigorous formulation of 
the assertion is the following: (l/^/2n)e~ x ' 2 Ax is approximately the probability for 



V n 

to belong to the interval (—Ax, Ax). 

Note that, in fact, p n is a random function, because it depends on the random variables 
X and S 2 . _ 

If n — > oo, then X — > x*, S 2 and /n — > 0, S 12 — > of yst in probability. Denote 



( X - X r 

Poo(x) = = e ^ 

T 2 

syst 



2na 2 



This is also a random function. The following proposition holds: 
Proposition. 

r b 

lim Pr{X (n) G [a, b}} = / Poo(a;)cix 

J a 

in probability, i.e., 

lim Pr (pr{X (n) G [a,b]} - [ p oc (a;)rfxl = 0, 
n ^°° I ia J 

if a = -(m — |), 6 = -(Z — |) /or some m, / G Z o£/ier words, a,b G -Z + 

This is a corollary of the limit theorems of probability theory (the law of large numbers 
and the central limit theorem) and condition (fT2l) . 

If we perform the repeated measurements, we can predict the probabilities of the results 
of the next measurement in the limit n — > oo using the constructed probability distribution 
function (fUj) . This justifies the use of the described construction. 



4 Conclusions 

It is shown that the use of real-valued probability density function as a fundamental concept 
of functional classical mechanics does not contradict to the thesis that the real irrational 
numbers are unobservable, since the density function is not a directly observable value. 
The construction of the probability density function based on the rational- valued results of 
measurements and an argumentation for this construction are given. It would be interesting 
to extend these results to the case of quantum mechanics, see [12]. 
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